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Abstract. 

The discovery of the first damped Lyman alpha (DLA) system in the early 1970s followed 
by the recognition that DLAs arise in intervening galaxies opened up a new field of galaxy 
evolution research. These highest HI column density absorption-line systems trace the bulk of 
the observed neutral gas in the universe, and therefore, have been used as powerful probes of 
galaxy formation and evolution back to the redshifts of the most distant quasars. The history 
and progress of DLA research over the past several decades is reviewed here. Larger datasets 
and deeper surveys, particularly over the last couple of years, have improved our knowledge of 
the neutral gas content and distribution in the universe at all observable redshifts, including 
the present epoch. New results on the statistics of DLAs at z < 1.65 from our HST-UV surveys 
are presented and discussed in the context of recent results at z — and at high redshift. We 
find that Qdla(z > 0) remains roughly constant to within the uncertainties; the z = value of 
the neutral gas mass density, Q g , is a factor of ~ 2 less than £Idla- The DLA incidence, n(z), 
undergoes rapid evolution between redshifts 5 and 2, but is consistent with the no-evolution 
curve in the current concordance cosmology for z < 2. We also show that if the local Schmidt 
law relating surface density of gas and star formation rate (SFR) is valid at the DLA redshifts, 
then the DLA SFR density is too low for them to provide a significant contribution to the cosmic 
star formation history (SFH) at z > 1. This implies that the DLAs are unlikely to be the same 
population as the star forming galaxies (i.e., the Lyman break and sub-millimeter galaxies) that 
dominate the SFH of the high redshift universe. We suggest that this discrepancy and the DLA 
"missing metals" problem could be the result of missing very high column density gas due to 
its very small absorption cross section. 
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1. Introduction 

The first published spectrum of a damped Lyman a (DLA) system, the z — 2.309 
absorber towards PHL 957 by Lowrance et al. (1971, 1972), was quickly followed by 
Beaver et al.'s (1972) independent observation of this object and their interpretation 
that the strong absorption line might be due to radiation damping from very high column 
density HI. They showed that the line profile was consistent with a column density of 
N(HI) = 2 x 10 21 atoms cm -2 , and noted its similarity to Lya absorption features along 
lines of sight to Galactic stars (Savage & Jenkins 1972). Serendipitous discoveries of DLAs 
towards the quasars Q1331+170, PKS 1157+014, and PKS 0528-250 were reported 
shortly thereafter (Carswell et al. 1975; Wright et al. 1979; Smith et al. 1979). The 
discovery of DLAs in the early 1970s in combination with the realization that Lya forest 
and narrow metal absorption lines arise in intervening material (Wagoner 1967; Bahcall 
& Spitzer 1969; Lynds 1971; Weymann et al. 1979; Sargent et al. 1980), led to a whole new 
approach to the study of galaxy evolution. By tracking absorption line properties with 
redshift, it was now becoming possible to explore the evolution of gaseous structures in 
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the universe. DLAs are the highest column density neutral hydrogen absorption lines seen 
in quasar spectra (N(HI) > 2 x 10 20 cm -2 ), and consequently, have been described as 
harboring the bulk of the neutral gas in the universe. Thus, barring any selection effects, 
tracking their properties as a function of redshift directly traces the cosmic evolution of 
neutral gas in the universe, and since cold neutral gas is a requirement for subsequent 
star formation, DLAs were soon recognized as important probes of galaxy formation and 
evolution (Wolfe, Turnshek, Smith, & Cohen 1986). As we describe below, Art Wolfe and 
collaborators were the first to exploit the potential of DLAs as cosmological probes, and 
in the early 1980s, began the first systematic search for DLAs. 

Prior to the deep imaging surveys of the last decade, DLAs were the only tracers of 
galaxies at high-redshift. A DLA system was also considered to be an unbiased probe since 
its detection was independent of the luminosity of the absorbing galaxy, and therefore, 
DLA-selected samples were thought to be more representative of the galaxy population. 
The detection of a DLA system only required the presence of a background quasar of 
sufficient brightness so that its spectrum could be obtained. Moreover, the strong line 
with damping wings characteristic of a DLA system could be detected in low-resolution 
spectra with signal-to-noise ratios not much greater than a few. In other words, they 
were ideally suited for searches through spectroscopic surveys of quasars. 

2. DLA Surveys: A brief history 

Since the redshifted Lya line can be observed with groundbased telescopes only for 
redshifts z > 1.65, initial DLA studies were restricted to this regime. Redshifted 21 
cm observations of radio loud quasars have also been used to detect large columns of 
neutral gas in absorption (see Briggs 1999 for an excellent review). Until the advent of 
UV telescopes, 21 cm absorbers along the line of sight to radio loud quasars were the only 
known high N(HI) systems below redshift 1.65. Discovery of the first redshifted 21 cm 
absorption line was reported by Brown & Roberts (1973) towards 3C 286 at z = 0.692, not 
long after the first optically discovered DLA, in a blind survey towards 100 (radio loud) 
quasars. The DLA towards 3C 196 was also discovered in this survey. Several more were 
detected towards quasars that had optically-identified low-ionization absorption lines of 
Mgll (Briggs & Wolfe 1983; see Briggs 1988 for a summary). By the early 1980s the pace 
of DLA discoveries was picking up, but it was clear that surveys were required for DLAs 
to become useful as statistical probes of neutral gas in the universe. In this section, we 
describe the progress of DLA studies beginning with optical surveys. It was also obvious 
that since our end point for galaxy evolution was all around us at z — 0, the neutral gas 
distribution at the present epoch could be used to anchor the evolutionary properties of 
galaxies observed in the distant universe. Rao & Briggs (1993) made the first attempt 
to describe the neutral gas distribution at z = in a framework for comparison with 
high redshift DLA studies. In §2.2, we describe this and subsequent z = surveys for 
"DLA" gas. Finally, in §2.3, we present results from surveys for z < 1.65 DLAs, the most 
recently explored redshift regime. 

2.1. Optical surveys for DLAs 

The Lick DLA Survey was the first systematic search for DLA candidates at redshift 
z ~ 2 (Wolfe et al. 1986). Wolfe et al. characterized the survey as a search for high 
redshift HI disks, motivated by the fact that DLAs had HI column densities comparable 
to those found in the Milky Way and in nearby spiral galaxies. The similarities to HI disks 
extended to their 21 cm properties, i.e., low spin temperatures, high 21 cm optical depths, 
and low velocity dispersions. The DLAs also showed low-ionization metal- line transitions 
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such as Mgll, Fell, CII, Sill, 01, etc., characteristic of cold clouds in the Milky Way's 
ISM. The survey included 68 quasars whose spectra were obtained at Lick Observatory at 
10 A resolution. This resulted in 47 absorption lines tagged as DLA candidates in a total 
redshift path of Az — 55. While this was an unexpectedly high incidence of DLAs given 
the then known incidence of Lya forest systems (Sargent et al. 1980), follow-up higher 
resolution spectroscopy (1-2 A) confirmed 15 bonafide DLAs (Turnshek et al. 1989; Wolfe 
et al. 1993). This led to the first determination of the incidence of DLA systems, i.e., the 
number of DLAs per unit redshift. They found dn/dz — 0.29 ± 0.07 at (z) = 2.5. 

This was followed by a larger survey of 100 additional quasars by Lanzetta et al. (1991). 
Including results from the Lick survey, they found 38 confirmed DLAs in 156 spectra over 
the redshift range 1.6 < z < 4.1 and a redshift path of Az = 155. With this larger DLA 
sample, they were able to determine for the first time the key statistical properties that 
are now routinely used to characterize the evolution of neutral gas, namely, the redshift 
number density dn(z)/dz [or n(z)], the cosmological neutral gas mass density £l g (z),\ 
and the column density distribution f(N) for the DLAs. Any theory of galaxy evolution 
could now be constrained by the observed properties of the neutral gas distribution in 
the universe. Lanzetta et al. found that the redshift number density could be expressed 
as a power law n(z) — rio(l + z) 1 with 7 = 0.3 ± 1.4, which implied that the incidence of 
DLAs did not evolve with redshift in the non-A cosmology that was assumed at the time. 
They also found that Q g (z = 2.5) « 1 x 10~ 3 was comparable to the mass density of stars 
in present-day spiral galaxies, and that it did not evolve over the redshifts probed by 
the survey. Their f(N) distribution was a power law, / (N) oc N~@, with j3 = 1.67, and 
showed an excess of systems above the DLA threshold in comparison to the extrapolated 
power-law distribution of lower N(HI) systems. 

Shortly thereafter, Wolfe et al. (1995) described results from spectroscopy of 228 
quasars drawn from the Large Bright Quasar Survey. Their statistical sample now in- 
cluded 58 confirmed DLAs at z > 1.65 in a redshift path length of over 300. In contrast 
to the Lanzetta et al. (1991) result, they found that Q g decreased with decreasing red- 
shift, and interpreted this evolution as due to depletion from star formation. In addition, 
they modelled the evolution in the f(N) distribution as being due to star formation in 
randomly oriented disks and found agreement with the observations if the SFRs followed 
the standard Schmidt-Kennicutt law seen in nearby galaxies. 

The next advancement in DLA studies from the ground occurred when surveys were 
extended to higher redshifts. As high redshift quasars became more numerous, so did high 
redshift DLAs; Storrie-Lombardi et al. (1996a, b), Storric-Lombardi & Wolfe (2000), and 
Peroux et al. (2001, 2003) extended this work beyond z = 4. Storric-Lombardi & Wolfe 
compiled a list of 81 DLAs in the redshift interval 1.65 < z < 4.69; they found that fl g 
was consistent with being constant for 2 < z < 4 and marginally consistent with a decline 
at z > 4. Peroux et al. discovered an additional 26 new DLAs, 15 at z > 3.5, in their 
survey of 66 z > 4 quasars. The highest redshift DLA now stood at z = 4.46. Peroux et 
al. also concluded that only 55% of the neutral gas mass at z > 3.5 was in DLAs and 
that the remaining was contained in Lyman limit and sub-DLA systems with HI column 
densities 10 19 < N(HI) < 2 x 10 20 cm~ 2 . Their result that the total neutral gas mass 
is conserved over redshifts 2<z<5 led them to suggest that the z « 3.5 turnover in Q, g 
was a signature of the epoch of formation and initial collapse of DLAs. 

With the advent of the Sloan Digital Sky Survey (SDSS) and its database of tens of 
thousands of quasar spectra, the next logical step was to mine it for DLAs. Prochaska 

f Q g = 1.3f2ff/, where the factor 1.3 corrects for a neutral gas composition of 75% H and 
25% He by mass. 
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& Herbert-Fort (2004) used the SDSS Data Release 1 to find 71 DLA systems in 1252 
quasar spectra at a mean redshift of z « 2.5. Contrary to the result of Peroux et al., 
they found that Lyman limit systems contribute less than 15% to f2 g at z > 3.5. This 
conclusion was based on the detection of 6 new DLAs at z > 3.5, half of which have 
N(HI) > 1 x 10 21 cm~ 2 . However, as noted by the authors themselves, the SDSS is most 
sensitive to DLA redshifts between 2.1 and 3, and the statistical uncertainties beyond 
z = 3.5 are still too high to make a definitive statement about the relative contributions 
of DLAs and sub-DLAs to the highest redshift values of Cl g . Prochaska & Herbert-Fort 
also found that n(z) is consistent with previous estimates, and that fl g decreases with 
decreasing redshift from z ~ 4 to z ~ 2 with a possible decline beyond z ~ 4. However, 
all of the data points in their figure 5 are within ~ la of each other, implying that 
fig is constant within this redshift range. As we show in §2.3, the low redshift data are 
consistent with this interpretation. Now, with their larger SDSS DR3 sample in which 
they detect over 700 DLAs beyond z — 2.1, Prochaska et al. (2005) report evidence for 
a statistically significant decline in f2 s from z ~ 4 to z ~ 2. The reader is referred to J. 
Prochaska's contribution in these proceedings for details on their new results. 

Fall & Pei (1989) first suggested that dust in high column density absorbers would 
influence the statistics of DLAs. The concern was that optically selected quasar samples 
would miss DLAs due to the reddening and dimming caused by foreground DLA host 
galaxies (see also Pei et al. 1999). To better evaluate this effect, Ellison et al. (2001) 
conducted a survey for DLAs towards radio-selected quasars and obtained spectra of 
all quasars regardless of their optical magnitude. They found that at the mean redshift 
of their sample, (z) = 2.4, both Udla and Qdla were only marginally higher than 
those measured using optically selected samples. Based on the uncertainties in their 
measurements, which were dominated by sample size, they concluded that £Idla could 
be no more than a factor of two higher than that derived from optical surveys, and that 
there was no hidden population of DLAs that was being missed in optical surveys. The 
effects of dust and other biases that might affect DLA statistics will be discussed further 
in §4. See Turnshck et al. (2005, these proceedings) for more on DLA selection effects. 



2.2. The HI distribution at z = 

The first attempt to describe the HI distribution at the present epoch in a form suitable 
for comparison with DLAs was made by Rao & Briggs (1993) and Rao (1994) who 
used HI 21 cm emission observations of a complete sample of nearby spiral galaxies 
to estimate f(N). The 27 galaxies used in the study included all galaxies with optical 
diameters greater than 7' that were in the declination range observable at Arecibo. They 
also excluded galaxies in the local group. This Arecibo survey (Briggs et al. 1980) was 
initiated to explore the distribution of HI in the outskirts of local galaxies with the 
primary aim of determining the extent of local galaxy halos. The results were then used 
to interpret quasar absorption line properties. For example, Briggs et al. argued that 
the lack of detectable HI beyond 2-3 Holmbcrg radii meant that local galaxies could 
not reproduce the observed high-redshift Mgll absorption line statistics. A decade later, 
we used these same data on the radial HI distribution of local galaxies to calculate 
f(N,z = 0) as described in Rao & Briggs (1993). In brief, if <A(M, N)> AN is the mean 
cross-sectional area subtended by gas with HI column density between N and N + AN 
in a galaxy of magnitude M, then 

c J2HM)(A(M,N))ANAM 
J[ ' H AN 
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Here, $(M) is the luminosity function of gas-rich galaxies and is used to normalize f(N). 
Since the optical luminosity of each galaxy in the sample was known and its radial N(HI) 
distribution was measured, (A{M, N)) AN, the area averaged over all inclinations, could 
be determined. A comparison with the (z) — 2.5 results of Lanzetta et al. (1991) showed 
that f(N,z = 0) was systematically lower than f(N,z = 2.5), and n(z = 0), which is the 
integral of f(N) over all DLA column densities, was a factor of «5 lower than n(z = 2.5). 
The cosmological mass density of HI, Qhi, which is proportional to / Nf(N)dN, was 
in agreement with other estimates of the z = mass density of neutral gas (e.g., Fall & 
Pei 1993), and was a factor of 3 less than the z = 2.5 value. Thus, we had found clear 
evidence for evolution of neutral gas properties to the present epoch. 

However, the local galaxy luminosity function as determined prior to 1993 did not 
include gas-rich galaxies that are now known to occupy the low-luminosity tail of the 
galaxy luminosity function. The overall normalization might have been underestimated 
as well. This offset manifests itself in the values derived for f(N), n(z), and Qhi- Ad- 
ditionally, the small sample and the large Arecibo beam size meant that improvements 
could be made (and were necessary) over this initial study. Moreover, the question of 
whether optically luminous galaxies contained all the HI was still being debated, and the 
possibility that HI clouds with no optical counterpart could contribute to the HI mass 
and cross section in the local universe was not completely discounted. To investigate this 
and to determine f(N) more accurately, many blind HI surveys were launched in the 
1990s. The Arecibo HI Strip Survey (AHISS; Sorar 1994; Zwaan ct al. 1997), the Arecibo 
Slice survey (Spitzak & Schneider 1998), the Arecibo Dual Beam Survey (Rosenberg & 
Schneider 2000), the Parkes all sky southern survey, HIPASS, (Zwaan et al. 2003, 2005; 
Ryan- Weber et al. 2003), and most recently, the Westerbork Synthesis Radio Telescope 
(WSRT) survey discussed by M. Zwaan in these proceedings have dramatically improved 
our knowledge of the HI distribution at the present epoch. Apart from determining f(N), 
n(z), and Hhi very accurately, these surveys showed convincingly that all of the neutral 
gas at the present epoch was accounted for by galaxies included in the optical luminosity 
function. Minchin et al. (2003), with HIDEEP - a 20 times more sensitive survey than 
HIPASS in a smaller region of sky - also found that all their HI detections had optical 
counterparts. This survey's limits, N(HI) « 4 x 10 18 cm~ 2 and M{HI) = 4 x 10 5 M Q , 
were the highest sensitivity limits achieved thus far. Thus, the consensus now is that 
there are no lurking HI giants or dwarfs in the local universe. While the HI mass density 
and cross-section are dominated by the high surface brightness, HI rich galaxies, Minchin 
et al. (2004) showed that low surface brightness galaxies (LSBs) contribute 30% to the 
HI mass density and 40% to the HI cross-section at the present epoch. This is very dif- 
ferent from what was assumed in the initial study of Rao & Briggs (1993), and as Zwaan 
(this proceedings) also asserts, is completely consistent with the high fraction of LSBs 
identified as DLA host galaxies at low rcdshift (e.g., Rao et al. 2003). 

2.3. The UV regime: z < 1.65 

Since the Lya line falls in the UV for redshifts z < 1.65, this most recent 70% of the age of 
the universe remained hidden until IUE and HST. The need for DLA surveys in the UV is 
illustrated in Figure 1 where £Idla is plotted as a function of time. (We now refer to f2 s 
as Qdla to specify the cosmological neutral gas mass density measured using DLAs. This 
change in notation is important because, as we discuss in §4, there is reason to believe 
that DLAs do not trace all the neutral gas in the universe.) As discussed above, both the 
optical, high-redshift and 21 cm, z = regimes have been serviced by all-sky surveys, 
namely, the SDSS and HIPASS. While an all-sky UV spectroscopic survey does not exist, 
and is unlikely, far-fetched, and wishful thinking on our part, we have, nevertheless, been 



6 



S. Rao 




I 0.5 



t/t 

Figure 1. The log of £Idla as a function of time (for the "737" cosmology: h — 0.7, Qm = 0.3, 
Qa = 0.7) with to representing z = 0. The open circle is Q g at z = from Zwaan et al. (2005) 
and the open square is the mass density of luminous matter at the present epoch from Panter 
et al. (2004). The high redshift (early times) data points are from Prochaska & Herbert-Fort 
(2004). The large gap in time accessible only by UV surveys for Lya illustrates the importance 
of the z < 1.65 redshift regime for DLA evolution studies. 

making incremental progress towards determining DLA statistics in this regime. The first 
concerted effort to obtain a large sample of UV spectra of quasars was made by the IUE 
satellite. Kinney et al. (1991) and Lanzetta, Turnshek, & Sandoval (1993) published a 
total of 260 IUE quasar and AGN spectra, which were subsequently used by Lanzetta, 
Wolfe, & Turnshek (1995) to estimate the incidence of DLAs in the UV. The IUE spectra 
proved to be unreliable, and follow-up spectra of the DLA candidates with HST confirmed 
only one of the four DLAs that were used to determine DLA statistics (Turnshek & Rao 
2002). Two were shown to not be DLAs, while the fourth was not observed with HST. 
Additionally, one system which was later shown to be a DLA did not satisfy the IUE 
sample's DLA criterion. 

DLAs at low redshift were turning out be very rare objects, and sample sizes had 
to be increased several fold to be able to get reasonably good statistics. Traditional 
blind surveys were not yielding useful results. In fact, the HST QSO absorption line Key 
Project found only one DLA in a redshift path of Az = 49 (Jannuzi et al. 1998). In 
Rao (1994) and Rao, Turnshek, & Briggs (1995), we used IUE and HST archival spectra 
to search for DLAs in previously-known, optically-identified Mgll systems. This was 
similar to the selection procedure used by Briggs & Wolfe (1983) to find 21 cm absorbers 
towards radio- loud quasars. Since low ionization metal lines are always seen along neutral 
gas sightlines in the Milky Way and in high-redshift DLAs, they act as tracers of high 
column density HI. Moreover, the Mgll doublet, AA2796,2803, is accessible for z > 0.11 
in optical quasar spectra. Beginning in the late 1980s, several groups conducted Mgll 
absorption line surveys and discovered many more Mgll systems than were available for 
the initial 21 cm searches of Briggs & Wolfe (Lanzetta, Turnshek, & Wolfe 1987; Tytler et 
al. 1987; Sargent et al. 1988, 1989; Caulet 1989; Steidel & Sargent 1992, SS92). The SS92 
survey, which was the largest, included 107 systems and resulted in the determination of 
the incidence of Mgll systems as a function of redshift and Mgll A2796 rest equivalent 
width. Thus, if a Mgll sample was selected appropriately and the fraction of DLAs in 
this sample could be determined, then the Mgll statistics could be used to calculate DLA 
statistics. We found 4 candidate DLAs in available IUE and HST archival data, and were 
able to place an upper limit on the values of n(z) and ^Idla at z = 0.8 by using this 
method. 
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In order to increase the low-redshift DLA sample, we obtained UV spectra of a sample 
of 60 additional Mgll systems with HST-FOS in Cycle 6. The sample of Mgll lines for 
our initial archival and Cycle 6 surveys was culled from the literature (Rao & Turnshck 
2000, RT00) and had a rest equivalent width limit of W^ 2796 > 0.3 A. The total sample 
of Mgll systems with UV Lya information now included 82 systems of which 12 were 
DLAs. With these results, we showed that Qdla was high at redshifts 0.1 < z < 1.65, 
consistent with a constant value for the neutral gas mass density from z w 4 to z w 0.5. 
This was primarily due to the high fraction of N(HI) > 10 21 cm~ 2 systems found in the 
survey. Our results also suggested that systems with both Mgll A2796 rest equivalent 
width, W^ 2796 > 0.5 A and Fell A2600 rest equivalent width, H 7 ^ 2600 > 0.5 A had a 
50% chance of being a DLA. Based on this result, we conducted a similar survey of 54 
Mgll systems in 37 quasars with HST-STIS in Cycle 9. Most of these satisfied the strong 
Mgll- Fell criterion for DLAs. Twenty seven had useful UV spectra and four of these were 
DLAs. The DLA towards Q1629+120 was discovered in this survey and was reported in 
Rao et al. (2003). The others are reported in Rao, Turnshek, & Nestor (2005a, RTN05). 
The strong Mgll- Fell criterion has since been modified using our new HST Cycle 11 
survey - see below. Meanwhile, using the same technique, Lane (2000) surveyed radio 
loud quasars for 21 cm absorption corresponding to 62 Mgll systems with the WSRT, 
and derived a similar result for SIdla at z < 1.65 (see Lane & Briggs 2001). 

With the number of quasar UV spectra increasing rapidly in the HST archives, it 
became possible to conduct a UV survey for Mgll at redshifts < z < 0.15. Churchill 
(2001) used 147 HST archival quasar/AGN spectra to obtain a total redshift search path 
of 18.8, and found 4 strong Mgll systems. For Wq 2796 > 0.6 A, he derived dn/dz = 
0.22lg Qg at (z) = 0.06, consistent with the extrapolation of dn/dz to lower redshifts 
for Wq 2796 > 0.6 A from Nestor (2004). All four quasars now have HST UV spectra 
that include redshifted Lya at the absorber redshifts. Three of the four have sub-DLA 
column densities and no information is available for the fourth due to an intervening 
Lyman limit system (Churchill 2001; Rao HST GO program 9382; Keeney et al. 2005). 
Thus, the number of DLAs in this sample is < 1. This gives a la upper limit to the 
DLA redshift number density of n(z = 0.06) ^ 0.19. While this result does not offer 
significant constraints, it is, nevertheless, consistent with the data points at z — and 
at low redshift (see §3). 

3. New results at low redshift 

Further progress could only be made with a dramatic increase in sample size. The 
SDSS sample of quasars, which numbered in the thousands when this phase of our Mgll- 
DLA project began, presented an unprecedented leap in the number of available survey 
quasars. The previous largest Mgll survey by Steidel and Sargent (SS92) used a sample 
of 103 quasars; the SDSS Early Data Release included nearly 4000. Nestor (2004) used 
SDSS-EDR quasar spectra to search for Mgll systems with the aim of quantifying the 
statistical properties of a large Mgll sample (Nestor, Turnshek, & Rao 2005, henceforth 
NTR05) and to conduct follow-up work to search for DLAs. In Cycle 11 (PID 9382), we 
targeted a sample of 83 Mgll systems with W A2796 > 1 A in 75 SDSS quasars. There 
were an additional 16 weaker Mgll systems observable in the same set of spectra. Overall, 
useful UV information was obtained for 88 systems, 25 of which are DLAs. For details 
of the sample and the selection criteria used to include Mgll systems in the survey, the 
reader is referred to RTN05. In total, we now have 197 z < 1.65 Mgll systems with 
Wq 2796 > 0.3 A for which UV observations of the Lya line have been obtained. Forty 
one of these are DLAs. 
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Figure 2. Plot of Wq 2796 vs. log N (HI) from RTN05. Filled circles are DLAs with 
N(HI) > 2 x 10 20 cm -2 . Arrows are upper limits in N(HI). Typical uncertainties are given by 
the error bars in the top left corner. 
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Figure 3. Distribution of Mgll A2796 rest equivalent widths, W^ 2796 . The gray histogram 
represents systems that are DLAs. Note that there are no DLAs in the first bin, i.e., for Mgll 
Wo 2796 < 0.6 A. The number of DLAs is roughly constant for Wq 2796 > 0.6 A, implying that the 
fraction of Mgll systems that are DLAs increases with increasing Wq 2796 . See RTN05. 



3.1. Metal line and HI correlations 

Since the systems in our sample were selected based on the rest equivalent width of Mgll 
A2796, Wo A2796 and N(HI) measurements exist for all 197 systems. Mgll A2803, the 
weaker member of the doublet, was also measured for all systems, while the Fell A2600 
and Mgl A2852 lines were measured for only a subset. Here, we explore correlations among 
metal line rest equivalent widths and HI column density. Figure 2 is a plot of Wq 2796 
versus log N (HI). We note that the upper left region of the figure is not populated 
implying that systems with Wq 279& > 2.0 A always have HI column densities N(HI) > 
1 x 10 19 atoms cm~ 2 . Figure 3 gives the distribution of Mgll W A2796 ; the DLAs form 
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Figure 4. Both histograms show the fraction of Mgll systems that are DLAs as a function of 
Mgll Wq 2796 , with the scale shown on the left axis. The solid circles are the logarithm of the 
mean HI column density in each bin for all observed Mgll systems in the left panel and only for 
the DLAs in the right panel. The scale is shown on the right. 



the gray histogram. It is noteworthy that there are no DLAs with Wq 2796 < 0.6 A.f 
In addition, the number of DLAs remains roughly constant for Wq 2796 > 0.6 A. As a 
consequence, the fraction of systems that are DLAs increases with increasing Wq 2796 . 
This is shown as a histogram in Figure 4; the y-axis on the left gives the fraction of DLAs 
as a function of Wq 2796 . We also plot the mean HI column density in each bin as solid 
circles with the scale shown on the right. Upper limits are assumed to be detections. 
The left panel includes all observed Mgll systems and the right panel includes only the 
DLAs. The vertical error bars are standard deviations in the mean and are due to the 
spread of N(HI) values in each bin, and the horizontal error bars indicate bin size. For 
the Mgll systems, there is a dramatic increase of a factor of ~ 36 in the mean HI column 
density from the first to the second bin, beyond which (N(HI)} remains constant within 
the errors. For systems with 0.3 A < W A2796 < 0.6 A, (N(HI)) = (9. 7 ±2.2) x 10 18 cm" 2 , 
and (N(HI)} = (3.5 ± 0.7) x 10 20 cm~ 2 for systems with W A2796 > 0.6 A. The right panel 
shows a trend for decreasing DLA column density with Wq 2796 . The reasons for this are 
not obvious, but arc likely to be due to small number statistics (see Figure 2), a real 
physical effect or a selection effect that is not yet understood (Turnshek et al. 2005). 

Figure 5a is a plot of Wq 2796 vs. Wq 2600 for systems with measured values of Wq 2600 , 
including upper limits. In RT00 we found that 50% of the 20 systems (excluding upper 
limits and 21 cm absorbers) with both W^ 2796 > 0.5 A and W^ 2600 > 0.5 A are DLAs. 
Now, with the expanded sample that includes 106 systems in this regime, we find that 36% 
are DLAs. The dashed line is a least-squares fit with slope b = 1.36. We note that DLAs do 
not populate the top left region of the diagram where the Wq 2796 to Wq 2600 ratio is > 2. 
In fact, if the sample is restricted to systems with Wq 2796 /Wq 2600 < 2, all but one of the 
DLAs from Figure 5a are retained, the outliers in the top left region are excluded as are 
most systems in the lower left corner of the plot. Figure 5b shows this truncated sample; 
the slope of the least-squares fit does not change significantly. We find b — 1.43. The only 
DLA that has been eliminated is the one with the smallest value of Wq 2600 . However, 

f Only one known DLA has lower metal-line rest equivalent widths. The 21 cm absorber at 
2 = 0.692 towards 3C 286 has Wo 2796 = 0.39 A and Wo A2600 = 0.22 A (Cohen et al. 1994). None 
of the 21 cm absorbers are included in our analysis here because they are biased systems (see 
RT00). 
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Figure 5 

of VK A260 °- Filled circles are DLAs. Typical error bars are 
dashed line is the best fit linear correlation with slope b 

Wo 2600 for systems with W^ 2796 /W^ 2600 < 2. Non-DLAs in the upper left and lower left regions 
of the diagram have been eliminated and the correlation is tighter than what is seen in Figure 
5a. Also, 38% of all systems are DLAs regardless of the values of Wo 2796 and Wo 2600 ; we believe 
this to be a more robust predictor of the presence of DLAs in Mgll-Fell systems. The slope of 
the best fit linear correlation is b — 1.43. Results are from RTN05. 



given the measurement errors for this system, its W^ 2796 /VF^ 2600 ratio is within la of 
2. This also implies that a system with metal line ratio W^ 2796 /W^ 2600 > 2 has nearly 
zero probability of being a DLA. For this truncated sample with W A2796 /W A2600 < 2, 
but no restrictions on the individual values of VF A2796 or VF A2600 , 38% are DLAs. In 
addition, all known 21 cm absorbers, including the z — 0.692 system towards 3C 286 
mentioned above, have VF A2796 /VK A2600 < 2. Thus, the W A2796 /W A2600 ratio is a more 
robust predictor of the presence of DLAs. 

This result is shown more dramatically in Figure 6. The ratio IF^ 2796 /!!^ 2600 is plotted 
as a function of N(HI) in Figure 6a. Ratios above 5 are not shown for clarity. These are 
mainly confined to log N (HI) < 19.2 with only one system above this column density at 
\ogN(HI) = 19.6. The DLAs populate the region of the plot where 1 <W A2796 /W^ 2600 < 
2; the two outliers lie within la of this range. A plot of the ratio W 7 ^ 2796 /!^ 2600 vs. Mgl 
^A2852 f Qr gygtems with measured values of IF^ 2852 , including upper limits, is shown in 
Figure 6b. Again, the DLAs are confined to the region where 1 <Wq 2796 /Wq 2600 < 2, but 
span the entire range of W A2852 . The two systems outside the range 1 < FF A2796 /tF A2600 < 
2 from Figure 6a do not have any information on Wq 28,52 . We find that 9 out of the 11 
systems with VF A2852 > 0.8 A are DLAs. We also note that systems with W A2796 /VF A2600 > 
2 are likely to have low values of W 7 ^ 2852 . 

For completeness, we also plot W A2600 vs. log N (HI) in Figure 7a and JF^ 2852 vs. 
log N (HI) in Figure 7b. There is no obvious trend in these distributions except for the 
fact that the upper left regions of the plots are not populated. There are no high rest 
equivalent width, low HI column density systems. This is not a selection effect since 
column densities as low as 10 18 cm~ 2 can often be easily measured. As is the case for 
W A2796 , this implies that systems with VF A2600 > lA or W> A2852 > 0.5A generally have HI 
column densities N(HI) > 10 19 cm~ 2 . Below this fairly sharp boundary, metal-line rest 
equivalent widths span all values of HI column density. 
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Figure 6. Left, a: Plot of Wq 2 '" 30 /Wq 2600 vs. N(HI). Systems with W^ 2796 /Wo 2600 > 5 are not 
shown for clarity; all of these have log N(HI) < 19.6. The DLAs are confined to the region of 
the plot where 1 <W X27m /W XW00 < 2. Right, b: Plot of W X27m /W ( X2600 vs. Mgl W X2852 . Solid 
circles are DLAs. Typical error bars are shown in the lower right corner. Again, the DLAs are 
confined to the region where 1 <W / A2796 /W A2600 < 2, but span the entire range of W A28 ° 2 . 
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Figure 7. Left, a: Plot of Wq 2600 vs. log N (HI). Arrows indicate upper limits. Typical uncer- 
tainties are given by the error bars in the top left corner. Right, b: Plot of Wo ° vs. log N(HI). 
Arrows indicate upper limits. Typical uncertainties are given by the error bars in the top left 
corner. 



3.2. Interpretation of Absorption Characteristics 

How can these trends be interpreted? Apart from the upper envelopes in Figures 2 and 7, 
there is no correlation between metal-line rest equivalent width and HI column density. 
Since the metal lines are saturated, the rest equivalent width is more a measure of 
kinematic velocity spread and not column density. High resolution observations of Mgll 
absorption lines have shown that the stronger systems break up into many components 
(e.g., Churchill et al. 2003), and span velocity intervals of up to 400 km/s. Turnshek et 
al. (2005) show line equivalent widths in velocity units of > 800 km/s in the strongest 
systems found in the SDSS. These highest equivalent width systems may arise in galaxy 
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groups; however, the more common systems like those in our DLA survey are more likely 
to arise in clouds that are bound in galaxy-sized potentials. A DLA is observed if one 
of the clouds along the sightline happens to be cold (less than a few 100 K), dense, and 
with a velocity dispersion of a few 10s of km/s. The more the number of clouds along 
the sightline, the higher the probability of encountering a DLA. This would explain 
the higher fraction of DLAs among large Wq 2796 systems and the lack of a correlation 
between Wq 2796 and N(HI). Only very rarely would a sightline intersect a single cloud 
resulting in small Wq 2796 and high N(HI), as in the 3C 286 system described in §3.1. This 
probabilistic approach to explain metal-line and HI strengths in high-N(HI) absorbers 
was also proposed by Briggs & Wolfe (1983) to explain their Mgll survey for 21 cm 
absorbers. They proposed a two-phase model where the 21 cm absorption is produced in 
galaxy disks, and the metal-line components that do not produce 21 cm absorption are 
produced in galactic halos. However, this multi-component/cloud model could be valid 
in any gas-rich galaxy, as is evidenced by DLA galaxy imaging studies (Le Brun et al. 
1997; Rao & Turnshek 1998; Turnshek et al. 2001; Rao et al. 2003; Turnshek et al. 2004). 
The disk models of Prochaska & Wolfe (1997) and the Haehnelt et al. (1998) models of 
infalling and merging clouds could reproduce these observations equally well. In other 
words, DLAs arise in pockets of cold gas embedded within warm diffuse gas or gas clouds 
in any bound system. 

Twenty one cm observations of low-redshift DLAs also reveal some cloud structure. 
For example, the z — 0.313 system towards PKS 1127—145 shows 5 components and 
the z — 0.394 system towards B0248+430 is resolved into 3 components (Lane 2000; 
Lane & Briggs 2001; Kanekar & Chengalur 2001). Since the Mgll line for these systems 
has not been observed at a resolution as high as the 21 cm observations, a one-to- 
one correspondence between the metal-line and 21 cm clouds cannot be drawn. In other 
instances, both warm and cold gas have been detected in a 21 cm DLA; Lane et al. (2000) 
find that two-thirds of the column density in the z = 0.0912 DLA towards B0738+313 
is contained in warm phase gas, and the rest is contained in two narrow components. 
The z = 0.2212 absorber towards the same quasar was also found to exhibit similar 
characteristics (Kanekar et al. 2001). In each of these cases, the line of sight probably 
intersects two cold clouds in addition to warm diffuse gas spread over a wider range of 
velocities that can be detected only in 21 cm observations of very high sensitivity. There 
are also several instances of DLAs not being detected at 21 cm (Kanekar & Chengalur 
2003). High spin temperatures (T s > 1000 K) corresponding to warm diffuse gas and/or 
covering factors less than unity towards extended quasar radio components have been 
suggested as possible explanations (Kanekar & Chengalur 2003; Curran et al. 2005). 

Clearly, a wide variety of cloud properties and their combinations are responsible for 
the observed properties of Mgll, DLA, and 21 cm absorption lines. Observations in the 
21 cm line have the added complication that the background radio and optical sources 
may sample different sightlines, thus making interpretation difficult. DLA observations, 
on the other hand, do not suffer from this effect. Large simulations of galaxy sightlines 
with varying cloud properties that reproduce the metal-line, DLA correlations shown 
in Figures 2-7 should be the next step towards improving our understanding of these 
absorption line systems. The simulations should not only reproduce the frequency of 
occurrence of DLAs in Mgll systems, but also the number density evolution of Mgll 
systems and DLAs. 
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3.3. The redshift number density of DLAs, ti^la 

The redshift number density of DLAs, udla, sometimes written as dn/dz, can be deter- 
mined using the equation 

riDLA(z) = f(z)n Mg ii(z), (3.1) 

where f(z) is the fraction of DLAs in a Mgll sample as a function of redshift and nMgii(z) 
is the redshift number density of Mgll systems. With the systems detected in SDSS EDR 
quasar spectra, Nestor (2004) derived 

dn/dz = N*(l + z) a e~^^ +z ^\ (3.2) 

where we have retained his notation and W = Wq 2796 (see also NTR05). N*, W* , a, and 
are constants. This expression is an integral over all Wq 2796 greater than Wq- Since 
our Mgll sample was assembled under various selection criteria that were modified based 
on our initial results (RT00), njjLA(z) had to be evaluated carefully. Again, we refer the 
reader to RTN05 for details on how Udla(z) was determined. 

Figure 8 shows the results for Udla{z) at low redshift split into two redshift bins (solid 
squares). We find 18 DLAs in 104 Mgll systems in the redshift interval 0.11 < z < 0.9 
with n DLA (z = 0.609) = 0.079 ± 0.019 and 23 DLAs in 94 Mgll systems in the redshift 
interval 0.9 < z < 1.65 with n DLA (z = 1.219) = 0.120 ± 0.025. The points are plotted at 
the mean redshift of the Mgll samples. The high- redshift data points are from Prochaska 
& Herbert-Fort (2004) and the z — point was estimated by Zwaan (in preparation) from 
a WSRT survey of HI in the local universe. The solid curve in Figure 8a is a no-evolution 
curve in the standard ACDM cosmology, that we refer to as the "737" cosmology, where 
(h,H,M = (0.7, 0.3, 0.7). This curve, which has been normalized at the z = data 
point, shows that the comoving cross-section for absorption declined rapidly by a factor 
of m 2 until z « 2 and has remained constant since then. This behavior might be a 
consequence of what has been observed in other studies of galaxy evolution, namely, that 
today's galaxies were in place by z w 1 and are a consequence of rapid merger and/or 
collapse events that occurred prior to this epoch. 

It has been customary in quasar absorption line studies to plot the logarithm of the 
redshift number density in order to illustrate its power law dependence with redshift, i.e., 
n(z) = n (l + z) 7 . In A = cosmologies, the exponent is a measure of the evolution in this 
quantity. For example, 7 = 1 for qo = or 7 = 0.5 for q = 0.5 imply no intrinsic evolution 
of the absorbers, and any significant departure from these values for 7 is evidence for 
evolution in the comoving number density of absorbers. We plot \ognDLA(z) as a function 
of log(l + z) in Figure 8b. The straight line is the power law fit to the data points 
with slope 7 = 1.27 ±0.11, and the curve is the same no-evolution function shown in 
the left panel. Thus, in the past the observations would have been interpreted as being 
consistent with the DLA absorbers undergoing no intrinsic evolution in a qo = universe, 
and marginally consistent with evolution in a qo — 0.5 universe. With the now widely 
accepted concordance cosmology, the interpretation has changed quite dramatically; as 
noted above, the nature of the evolution is redshift dependent. 

The possibility of a dust bias in DLA statistics at low redshift was addressed by Ellison 
et al. (2004) using the same CORALS survey that Ellison et al. (2001) used to show that 
the effect of dust at z m 2 was minor. They determined the Mgll dn/dz for redshifts 
0.6 < z < 1.7, and found excellent agreement with the results of NTR05 who used the 
SDSS-EDR quasar sample. Based on the RT00 strong Mgll-Fell selection for DLAs, they 
concluded that dust bias was not important for DLAs at low redshift either. 
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Figure 8. Left, a: Plot of Udla(z) versus redshift. The new low-redshift data points from 
RTN05 are shown as filled squares. The high-redshift points (filled circles) are from Prochaska 

6 Herbert-Fort (2004) and the z = data point is from an analysis of local HI (Zwaan, in 
preparation). The solid line is a no-evolution curve in the "737" cosmology with (h,Q.M ,&a) = 
(0.7, 0.3, 0.7) normalized at z = 0. The comoving cross-section for absorption declined rapidly 
by a factor of m 2 until z fs 2 and has remained constant since then. Right, b: Plot of log udla{z) 
as a function of log(l + z). The straight line is the power law fit to the data points with slope 

7 = 1.27 ± 0.11, and the curve is the no-evolution function shown in Figure 8a. 



n DLA (z) = ^^L nDLA {z) (N(HI)) 7^^, (3.3) 



3.4. The cosmological mass density Qdla 
We can determine firjLA from the DLA column densities and tidi,a{z) via the expression 

pm H H /ATtTTTW E ( z ) 

n DL A{z) (N{HI)) 

cp c (1 + zy 

where 

E(z) = (fi M (l + zf + (1 - Q M - fi A )(l + zf + O0 1/2 - (3.4) 
Again, the "737" cosmology has been used in the calculation of Qdla- Also, [i = 1.3 
corrects for a neutral gas composition of 75% H and 25% He by mass, wih is the mass 
of the hydrogen atom, p c is the critical mass density of the universe, and (N(HI)) is the 
mean HI column density of DLAs in each redshift bin. 

In contrast to the redshift number density evolution shown in Figure 8, we find that 
la has remained constant from z = b to z = 0.5 to within the uncertainties. Figure 9 
shows the new results as solid squares. Specifically, for the redshift range 0.11 < z < 0.90, 
we find n DLA (z = 0.609) = (9.7 ± 3.6) x 10~ 4 and for the range 0.90 < z < 1.65 we get 
^dla( z = 1.219) = (9.4±2.8) x 10 -4 . The uncertainties have been reduced considerably in 
comparison to our results in RT00. The reasons for this are two fold. First, the uncertainty 
in riMgii has been significantly reduced due to the fact that the Mgll sample size was 
increased 10-fold. Second, the number of DLAs in each bin has increased by more than 
a factor of 3. Thus, the uncertainties in the low- and high-redshift data points are now 
comparable. Note that the statistics of the high-redshift data are also improved due to 
the inclusion of an SDSS DLA sample (Prochaska & Herbert-Fort 2004). Nevertheless, 
our basic conclusion from RT00 has remained unchanged, namely, that the cosmological 
mass density of neutral gas remains roughly constant from z ps 5 to z i=s 0.5. 

The drop in redshift number density from z — 5 to z = 2 along with a constant mass 
density in this range indicates that while the product of galaxy cross-section and co- 
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Figure 9. Cosmological mass density of neutral gas, Qdla, as a function of redshift. The filled 
squares are the new low-redshift data points from RTN05. The high-redshift points (filled circles) 
are from Prochaska & Herbert-Fort (2004) and the open circle at z = is from Zwaan et al. 
(2005). The open square at z = is the mass density in stars estimated by Panter et al. (2004) 
from SDSS data. While the statistics have improved considerably, our basic conclusion from 
RT00 has remained unchanged, namely, that the cosmological mass density of neutral gas has 
remained constant from z ~ 5 to z ~ 0.5. 



moving number density is declining, the mean column density per absorber is increasing. 
This is, again, consistent with the assembly of higher density clouds as galaxy formation 
proceeds. 

On the other hand, a constant cross-section from z k, 1 to z = along with a drop 
in mass density from z « 0.5 to z = is indicative of star formation that depletes the 
highest column density gas while keeping the absorption cross-section constant. This 
would in turn require that the column density distribution of DLAs change such that the 
ratio of high to low column densities decreases from low-redshift to z = 0. As we will see 
in the next section, the column density distribution does show some evidence for this. 

3.5. The column density distribution f(N) 

Figure 10 shows the normalized cumulative column density distribution (CDD) for the 
three redshift regimes. The dotted curve is the 2 = distribution derived by Ryan- Weber 
et al. (2003, 2005) from HIPASS data. The red (thick, solid) curve is derived from the 
41 DLAs at low redshift (RTN05) and the green (thin, solid) curve is from the "total" 
sample of Prochaska & Herbert-Fort (2004) . The change in the three CDDs with redshift 
is exactly what is expected based on the udla and £Idla results. Namely, that the 
low-redshift CDD shows a higher incidence of high column density systems than at high 
redshift presumably due to the assembly of gas as galaxy formation proceeds, followed 
by a decrease in the fraction of high column density systems to z = 0, presumably due to 
the depletion of gas during star formation. Thus, at least qualitatively, the evolutionary 
behavior of udla, ^dla, and the CDD are entirely consistent with one another. A KS 
test shows that there is a 25% probability that the high- and low-redshift curves are 
drawn from the same population; this is significantly higher than what we observed in 
RT00, where the two samples had only a 2.8% probability of being drawn from the same 
population. However, the general trend that the low-redshift sample has a higher fraction 
of high column density system still remains. 

Figure 11 is a plot of the log of the absolute CDD function, log /(TV), as a function of 
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Figure 10. The normalized cumulative CDD of DLAs for three redshift regimes. The red, thick 
solid line includes the 41 low-redshift DLAs at (z) = 0.95 (RTN05). The green, thin solid line 
includes 163 high-redshift DLAs at (z) = 2.94 (Prochaska & Herbert-Fort 2004), and the blue, 
dotted curve is the z — estimate from Ryan- Weber et al. (2003, 2005). 
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Figure 11. Absolute CDD function for the three redshift regimes. The red, thick solid line is a 
least squares fit to the low-redshift data points (RTN05) and has slope = 1.4; the green, thin 
solid line is a least squares fit to the high-redshift data points (Prochaska & Herbert-Fort 2004) 
with p = 1.8. The dotted line is the CDD with /3 = 1.4 for log N (HI) < 20.9 and (3 = 2.1 for 
log N (HI) ^ 20.9 (Ryan- Weber et al. 2003, 2005). 

log N (HI) for the three redshift regimes. The turnover with redshift is most apparent in 
the lowest and highest column density bins. We derive (3 = 1.4 ± 0.2 and (3 = 1.8 ± 0.1 
at low and high redshift respectively, where the CDD is expressed as f(N) = BN~@ '. 
At z = 0, Ryan- Weber et al. (2003) derive f3 = 1.4 ± 0.2 for log N (HI) < 20.9 and 
/3 = 2.1 ±0.9 for log N(HI) ^ 20.9. The general form of the absolute CDD does not vary 
considerably with redshift, which in turn explains the roughly constant value of £Idla- 
The differences in the f(N) distributions are subtle, implying that the gas content in 
DLAs is not changing drastically. This is strong evidence that DLAs do not have high 
SFRs and are, therefore, a different population of objects than those responsible for 
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much of the observed luminosity in the high rcdshift universe. On the other hand, a non- 
evolving DLA population might be observed if the gas that is used up in star formation is 
replenished from the inter-galactic medium at a comparable rate. This possibility seems 
rather contrived, and requires more proof than the current observational evidence can 
provide. In the next section, we discuss further evidence that suggests that DLAs and 
Lyman break galaxies (LBGs) are mutually exclusive galaxy populations. 



4. The Nature of DLAs 

With the identification of DLAs in quasar spectra, follow-up work involving the study 
of DLA galaxies, their environments, DLA gas metallicities, kinematics, temperatures, 
and ionization conditions, and numerical and semi- analytic modelling has kept many 
astronomers busy for decades. A keyword search for "damped Ly" in NASA's abstract 
database returns nearly 1000 results, too many to reference here. Despite this, a consensus 
on the nature of these objects is still lacking. One of the primary reasons for this is that 
the selection effects inherent in the observations are still being sorted out. Many of these 
issues were raised at this conference, and we address some of them here. See also Turnshek 
et al. (2005, these proceedings) for further discussion on DLA selection effects. 

The conclusion that DLA surveys identify the bulk of the neutral gas in the universe 
is based on three results or assumptions. First, integration of the CDD shows that a 
relatively small fraction of the neutral gas is contributed by Lyman limit and sub-DLA 
absorption systems with 3 x 10 17 < N(HI) < 2 x 10 20 cm~ 2 , at least for z < 3.5 (Peroux 
et al. 2003), and perhaps at all redshifts (Prochaska & Herbert-Fort 2004). Second, as 
discussed earlier, dusty regions do not cause DLA surveys to miss a large fraction of the 
neutral gas (Ellison et al. 2001, 2004). Third, the biases introduced by gas cross section 
selection are small. The important points to emphasize here are that the interception (or 
discovery) probability is the product of gas cross section times comoving absorber number 
density, and no DLAs with N(HI) > 8 x 10 21 cm~ 2 have been discovered. Thus, the third 
assumption requires that rare systems with relatively low gas cross section and very high 
HI column density arc cither absent or have not been missed to the extent that the neutral 
gas mass density will be significantly underestimated by quasar absorption line surveys. 
This third assumption might have to be reevaluated in order to explain the discrepancy 
between the star formation history (SFH) of DLAs as inferred from their HI column 
densities and that determined from galaxies that trace the optical luminosity function. 
We explored this connection in Hopkins, Rao, & Turnshek (2005) and summarize it here. 

We can use the global Schmidt law relating gas surface density and SFR to infer the 
star forming properties of the DLAs. Assuming that the Schmidt law, quantified by 
Kennicutt (1998), is valid at the DLA redshifts we can convert DLA column densities 
into gas surface densities and estimate their SFR. In Hopkins et al. (2005) we show that 
the SFR density of DLAs relates to udla{z) and VIdla{z) as follows: 



where p* is the SFR density in units of M yr _1 Mpc -3 and dX/dz is the absorption 
distance which equals (c/Hq)(1 + z) 2 /E(z). E(z), /x, and p c arc as defined in §3.4. Fig- 
ure 12 shows the SFR density of DLAs estimated using equation 4.1 and the values of 
nDLA{z) and £Idla(z) plotted in Figures 8 and 9. The hatched region encompasses the 
majority of SFR estimates in the literature (Hopkins 2004) and the solid curve is a fit 
to these data. The crosses are estimates of the SFR in DLAs using the CII* absorption 
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Figure 12. SFR density of DLAs (filled circles) as a function of redshift from Hopkins, Rao, 
& Turnshek (2005). Crosses are estimates from Wolfe et al. (2003) using the CII* method (but 
see footnote). The hatched region encompasses the majority of SFR estimates in the literature 
(Hopkins 2004) and the solid curve is a fit to these data. 



line technique of Wolfe et al. (2003). The CII* results are consistent with the fitted curve 
derived from the SFH of luminous galaxies, and have been interpreted as evidence that 
the DLA population is the same as the LBGsf. However, the few high-rcdshift DLAs 
for which host galaxies have been identified yield SFRs that arc at the very faint end 
of values determined for LBGs (e.g., Weatherley et al. 2005, M0ller et al. 2002). Recent 
results from semi-analytic models and numerical simulations by Okoshi & Nagashima 
(2005) and Nagamine et al. (2004) indicate that the average masses of DLA galaxies 
are likely to be low, around 10 9 M Q , which is consistent with low integrated levels of 
SFR. Okoshi & Nagashima (2005) find a broad range of SFRs from 10~ 6 to 100 M Q 
yr _1 , with a mean value of 0.01 M Q yr _1 for DLAs at redshifts z < 1, and conclude 
that DLAs are typically dwarf systems with low SFRs. Therefore, it is entirely possible 
that the high redshift observations of Weatherley et al. (2005), and M0ller et al. (2002) 
are sampling the brightest end of DLA galaxy properties. Lanzetta et al. (2002) have 
also shown, through similar HI surface density and SFR arguments that LBGs have SFR 
surface densities that are up to 3 orders of magnitude higher than that of the DLAs. In 
fact, their SFR densities do not overlap for the observed DLA redshift range 1.65 < z < 5 
(see Lanzetta et al. figure 3b). This implies that the high SFR surface density objects 
have very high column density neutral plus molecular gas (10 22 to 10 25 cm -2 ). These are 
comparable to the surface gas and SFR densities seen in local starburst nuclei (Kenni- 
cutt 1998). Observations of local star forming galaxies have shown that N(HI) is always 
> 10- 2 N(H 2 ) in the high SFR regions (e.g., Wong & Blitz 2002), and therefore, could lie 
above the currently observed DLA regime. If these systems have very small interception 
cross sections, they may be missed in surveys for DLAs. For example, an absorber with a 
size of about 100 pc, comparable to giant molecular clouds (GMCs) which are the sites of 
star formation, has a cross section that is « 10 4 times smaller than known DLAs, which 
typically have effective radii of « 10 kpc (Monier et al. 2005). Assuming that there are on 
the order of 10 GMCs per galaxy, the total cross section per unit volume, i.e., intercep- 
tion probability, for these very high column density gas systems would be on the order 

f Using more recent UV background radiation estimates, Wolfe (2005, these proceedings) 
finds lower CII* estimates of DLA SFR density which are now in agreement with the Hopkins 
et al. result. 
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of 10 3 times smaller. This means that 10 3 DLAs need to be detected in order to find one 
very high column density system. With the SDSS, we are getting close, but are not quite 
there yet. A one in a thousand system with N(HI + H2) = 10 24 cm~ 2 would increase 
the SFR density of DLAs by more than a factor of 2, and bring the DLA SFR density 
into agreement with the luminous SFR density. Searches for molecular gas in DLAs have 
resulted in only a handful of detections. Moreover, the molecular gas fraction in the few 
DLAs with H 2 detections is very small (e.g. Ledoux et al. 2003), and is consistent with 
the idea that the known sample of DLAs does not trace the majority of the star forming 
gas in the universe. 

A similar conclusion is reached when considering the low metallicities of DLAs. Figure 
13a shows the most recent compilation of DLA metallicities from Rao et al. (2005b). 
Pcttini and collaborators were the first to make metallicity measurements of a large 
number of DLAs; the early data indicated that DLA metallicities do not evolve from 
redshifts z w 3.5 to z « 0.5 (Pettini et al. 1999 and references therein). With a much larger 
sample, Prochaska et al. (2003) and Rao et al. (2005b), the latter with more measurements 
at z < 1.65, showed that DLA metallicities do show evidence for evolution. In Figure 
13b, we show the metal mass density as a function of redshift. The DLA metallicity data 
from Rao et al. (2005b) are converted to metal mass by assuming logS!^^ = —3 and a 
solar metal mass fraction of Zq = 0.02. The metal mass density in stars, shown as the 
hatched region in Figure 13b, can be derived from the SFR, since the SFR is related to 
the metal production rate (e.g., Madau et al. 1996). Stellar population synthesis results 
indicate that /?* = 63.7 pz (Bruzual & Chariot 2003, Conti et al. 2003). The z — and 
z = 2.5 results of Dunne et al. (2003) are also shown. Dunne et al. used the dusty, 
high SFR, submillimeter galaxy population at z = 2.5 to derive their result, which is 
reasonably consistent with the evolution of pz from the SFH. They concluded that high 
redshift DLAs are not the same population as these star forming galaxies, which contain 
the majority of the metals at this epoch. The DLAs have about two orders of magnitude 
less mass in metals than star forming galaxies. The combination of low metallicity and 
low gas mass density naturally leads to a space density of metals that is significantly 
lower than in the luminous galaxy population. This "missing metals" problem cannot be 
solved simply by invoking spheroids, dwarfs, low surface brightness galaxies, or gas at 
large galactocentric distance. Instead, this problem may in fact be the result of missing 
a substantial fraction of the metal-enriched gas in DLA surveys, for example in the 
GMCs mentioned above, due to their very small interception probability. That molecules 
are detected in DLAs with higher than average metallicities is consistent with this idea 
(Ledoux et al. 2003). Thus, the SFH and missing metals problem can both be attributed 
to this same selection effect. All of this holds true at high redshift. We note, however, that 
the SFR in DLAs and luminous objects at z < 1 are comparable, but their metal mass 
densities are not. This can be explained by noting that by z « 1, the SFR in luminous, 
high SFR, galaxies begins to decline, i.e., the luminosity density of the universe fades, 
leaving behind the low SFR DLA galaxies as the main contributor to the volume averaged 
SFR. The metal mass density, on the other hand, has been building up in the high SFR 
galaxies at high redshift, and continues to remain high. In fact, pz plateaus at z < 1 
because the DLAs have low SFRs and do not increase the metal mass density significantly. 

However, the possibility that these very high column density gas systems are being 
missed by DLA surveys may not only be due to their small gas cross sections, but 
also because they are likely to be very dusty. Ledoux et al. (2003) find that the DLAs 
in which H 2 is detected have among the highest metallicities and the highest depiction 
factors, hinting at the possibility of much higher depletions in much higher column density 
molecular gas clouds. These are not found in the radio loud quasar surveys for DLAs 
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Figure 13. Left, a: Metallicity vs. redshift for DLAs plotted as open squares where the area of 
each square is proportional to N(HI) (Rao et al. 2005b, Prochaska et al. 2003). The cosmic mean 
metallicities weighted by N(HI), {Z), are filled stars and the unweighted mean metallicities are 
filled circles. The best fit linear slope is m = —0.26 ±0.06, indicating increasing metallicity with 
time. Right, b: Metal mass density as a function of redshift from Hopkins et al. (2005). Solid 
circles are DLA metallicities from Figure 13a converted to mass density, and the hatched region 
is derived from the SFH corresponding to the hatched region in Figure 12. Triangles at z = 
are from Calura & Matteucci (2004); open circles are from Dunne et al. (2003). 



because these surveys also suffer from the small cross section selection effect. Not enough 
radio loud quasars have yet been surveyed to find the putative one in a thousand very 
high column density system. 

Gravitational lensing has the opposite effect on DLA surveys. Magnification by DLA 
galaxies could brighten background quasars, and preferentially include them in magnitude- 
limited samples. Lc Brun et al. (2000), with HST imaging observations, showed no ev- 
idence for multiple images of background quasars and concluded that the quasars were 
magnified by at most 0.3 magnitudes. In addition, Ellison et al. (2004) and Peroux et al. 
(2004) using statistical tests on low redshift Mgll and DLA samples, showed that lens- 
ing bias is a minor effect. More recently, using the SDSS Mgll survey results of Nestor 
(2004), Menard et al. (2005) show that quasars behind strong Mgll absorbers, of which 
DLAs are a subset, show little magnification bias, and that its effect on Qdla at low 
redshift is negligible (see also Menard 2005). It is also unlikely that the lowest redshift 
points that we derived from our HST-UV data (Figure 9) are affected by lensing bias. 
This is because the DLAs with the highest HI column densities at z rs 0.5 arise in dwarf 
galaxies (Rao et al. 2003), and consequently, do not have the mass required to produce 
significant magnification. 



5. Summary 

Since the discovery of DLAs in the early 1970s, a considerable amount of effort has 
been put into deciphering the properties of neutral gas in the universe. The observational 
techniques that have to be used to trace the neutral gas separate the Lya universe into 
three redshift regimes: the optical/IR at z > 1.65, the UV at < z < 1.65, and the radio 
(21 cm) at z — 0. The current status of DLA surveys in these three redshift regimes was 
reviewed; while the z = and z > 1.65 regimes have been more than adequately serviced 
by nearly all-sky surveys, the UV regime could do with an order of magnitude more 
DLAs. We also put forward a change in the paradigm that DLAs trace the bulk of the 
neutral gas in the universe. Selection effects that preclude us from observing the highest 
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column density gas might be affecting our interpretation of what DLA statistics really 
mean. Evidence for this comes from a comparison of DLA and luminous mass densities, 
SFRs, and metal mass densities as a function of redshift. See Turnshek et al. (2005, these 
proceedings) for more discussion on DLA selection effects. 
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